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1 Introduction 


Observational evidence from cosmology has confronted theorists with three 
crucial questions: 

• What drove primordial inflation? 

• What is causing the current phase of acceleration? 

• How can we dehne cosmological observables which are gauge indepen¬ 
dent, infrared hnite, and BPHZ renormalizable (Bogoliubov, Parasiuk, 
Hepp and Zimmerman [H El El 0]) in the sense of low energy effective 
held theory [5l [6]? 

This paper concerns the third question, whose answer in hat space would 
be either the S-matrix (when only massive particles are present) or else in¬ 
clusive rates and cross sections (when massless particles occur). Although 
a formal S-matrix can be constructed for massive scalars on de Sitter back¬ 
ground [7], the causal structure of this geometry precludes local observers 
from measuring this quantity, so it cannot serve as an observable. To the 
tree order accuracy which has so far been resolved [8], the primordial scalar 
power spectrum seems to be well-represented by a 2-point quantum 

gravitational correlation function. This correlator can be given a gauge in¬ 
dependent expression at tree order |9l [10], but no local extension of it can 
be gauge independent at higher orders [m. Hence dependence upon the 
gravitational gauge has emerged as a central issue in loop corrections to cos¬ 
mological observables [l2l[l3l[Il[l5l[l6l[I71[l8l[l9l[20ll^|22l[23l[Ml[25l^ 

[27l[28l[29l[3ni|3llE2lE3lEaE5]. 



Figure 1: Feynman diagrams relevant to the one loop vacuum polarization 
from gravitons. Wavy lines are photons and curly lines are gravitons. 

The aim of this paper is to explore gauge dependence in the one graviton 
loop correction to the vacuum polarization x') on de Sitter back¬ 

ground. This quantity, whose diagrammatic depiction is given in Fig. [H 
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explores the same issues of gauge dependence as but is quite a bit 

simpler to compute and renormalize. It can be used to quantum-correct 
Maxwell’s equation |36] . 


dy 




+ / (fx' 




{X]x')Aiy{x') = J^{x) 


( 1 ) 


where A^{x) is the vector potential, = d^A^, — dyA^ and J^{x) is the 
current density. Quantum corrections to dynamical photons emerge from 
solutions with J^{x) = 0, whereas quantum corrections to electrodynamic 
forces derive from the response to nonzero J^{x). 

We dehne the graviton held operator h^y{x) by subtracting the de Sitter 
background g^y{x) from the full metric held operator 

= gi,y{x) + Kh^y{x) , = levrG . (2) 


We shall compute the vacuum polarization in the 1-parameter family of exact 
covariant gauges that is a generalization of the de Bonder condition, 


-na 


^Dphpfj 


= 0 


6 > 2 , 


(3) 


where Dp stands for the covariant derivative operator in de Sitter back¬ 
ground, treating hpy as a 2nd rank tensor. In gauges of the sort (IHl) the 
graviton propagator breaks up into a transverse-traceless “spin two” park 
which does not depend upon 6, and a 6-dependent, “spin zero” part l38] J^l 


afi^'yS 


x: x') = i 


A2 

aP‘^'yS 


{x-,x)+i 


AO 


[X] X 


(4) 


The graviton contribution to the vacuum polarization on de Sitter must de¬ 
pend upon h because its hat space limit takes the form |3S], 




flat 


^^r(f)r(f-i) 

16(B —l)7r'^ 


C2 + C'o(^) 




1 

Ax'^^~‘^ 


(5) 


Because covariant gauge propagators derive from both constrained and dynamical 
fields, the spin two part of (jd]) — whose zero components do not vanish — actually comes 
from both the ^{D — 3)D dynamical gravitons and Z? — 1 of the constrained fields. The 
spin zero part derives from the remaining constrained field. 
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Here r]^‘' is the Lorentz metric and Ax"^ = r]fj^^{x — x')^{x — x'y. The gauge 
independent, spin two term contributes a coefficient which vanishes in D = A 
spacetime dimensions, 


6*2 


{D-A){D-2f{D + l){D + 2) 

4(£>-l) 


( 6 ) 


and the gauge dependent, spin zero coefficient is. 


Note that C'o(&) is negative semi-dehnite for D = A. 

It is commonplace to dismiss as unphysical any gauge dependent quan¬ 
tity like ([5]) |T6]; however, that view is simplistic. The gauge-independent 
S-matrix of flat space arises from combining gauge-dependent Green’s func¬ 
tions, so the latter must possess legitimate physical information mixed in with 
artefacts of gauge hxing nn. From the manner in which the flat S-matrix is 
constructed HDl, one realizes that this physical information is distinguished 
by possessing momentum space poles on each external leg. In position space 
these poles correspond to secular growth when the Green’s function is inte¬ 
grated against tree order mode functions, as [^n‘'](x; x') necessarily is in the 
perturbative solution of ([T]) for dynamical photons. It is therefore reasonable 
to expect that the leading secular growth factors of de Sitter Green’s func¬ 
tions might be gauge independent [2^. That is what we seek to check for 
[^n^](x; x'). We will of course be able to compare our results for different 
values of 6 > 2. We can also compare with the result previously obtained 
[H] in a noncovariant gauge |l2l 03] . 

Some comments on the physics are worthwhile before commencing this 
difficult computation. One might think an uncharged held like h^^{x) is not 
capable of contributing to vacuum polarization but this ignores the role of 
electric and magnetic helds in transferring momentum. Virtual gravitons 
which interact with photons — either real or virtual ones — can alter this 
momentum. There is no change in how dynamical photons propagate on hat 
space background [3H|, essentially because virtual gravitons ahect a single 
photon the same way throughout space and time. However, the interaction 
between charged particles on hat space background is slightly strengthened 
at short distances, as can be inferred from the gauge independent scattering 
amplitude [H]. One way to understand this ehect is that the virtual pho¬ 
tons which transfer momentum between nearby charges do not survive long 


Coih) = ~{D-2f 


fDb—2^ 


(Db-2\ 

2{D-AY 

1 6-2 J 

\D-2) 

V 6-2 J 
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enough to experience the full effect of buffeting by the longest wave length 
(and hence longest lived) virtual gravitons. 

On de Sitter background the vacuum polarization must of course show 
the same strengthening of force at short distance that is encoded in its flat 
space limit [39l 01]. However, it should also manifest new, secular effects 
arising from the inflationary production of gravitons. For example, a tree 
order photon redshifts as it propagates, whereas the continual replenishment 
of Hubble-scale gravitons should lead to a relative one loop enhancement, 
as momentum tends to flow from inflationary gravitons into the ever weaker 
photon. Similarly, the force between widely separated sources should be rel¬ 
atively enhanced because the highly infrared virtual photons which mediate 
the force are more likely to acquire momentum from, rather than lose it to, 
the constant pool of Hubble-scale gravitons. Both effects have been seen 
[I5l06] when the noncovariant gauge [^n^](x; x') [H] is used in equation ([1]), 
and it will be fascinating to learn what happens in our covariant gauge ([3|). 

The one-loop effects of inflationary gravitons in the noncovariant gauge 
021 03] have been studied for a variety of other particles over the years: 

• The graviton self-energy has been computed m, but has not yet been 
used to quantum-correct the linearized Einstein equation. However, the 
Hartree approximation has been used to show that the Weyl curvature 
of dynamical gravitons experiences a secular enhancement |48j . 

• The self-energy has been computed for massless fermions [IHI and, to 
hrst order in the mass, for massive fermions [50]. Quantum correcting 
the Dirac equation reveals a secular enhancement of the held strength 
of massless fermions [5I1152]. The result for the massive case has not 
yet been derived. 

• The self-mass of massless, minimally coupled scalars has been computed 
[53] . However, quantum-correcting the Klein-Gordon equation shows 
no secular enhancement of the scalar held strength [53]. 

• A partial result has recently been obtained for the self-mass of massless, 
conformally coupled scalars [55], but it has so far not been used to 
quantum-correct the linearized held equation. 

The fact that inhationary gravitons give secular enhancements to the held 
strengths of massless fermions, gravitons and photons, but not to massless. 
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minimally coupled scalars, seems to be due to spin [SS]. In each case energy 
and momentum of the physical particle under study redshifts as it propa¬ 
gates. If inflationary gravitons can only interact with the particle through 
its redshifting energy-momentum then the interaction cuts off rapidly and 
there can be no growing effect at late times. The presence of spin gives rise 
to a new interaction which does not cut off, so that particles with spin are 
scattered more and more as they propagate further through the sea of infla¬ 
tionary gravitons. Any other interaction which persists to late times should 
give rise to the same sort of secular enhancement, which is why it will be fas¬ 
cinating to see what happens to slightly massive fermions and to conformally 
coupled scalars. 

This paper consists of hve sections, of which the hrst is this Introduction. 
In section [2] we give those reductions of the diagrams in Fig. [1] which do 
not depend upon the form of the graviton propagator. Section [3] derives the 
contribution from the spin two part of the graviton propagator, and section 
m computes the contribution from the spin zero part. Our conclusions are 
presented in section |3 

2 Preliminary Reductions 

The purpose of this section is to describe those parts of the computation 
which do not require a specihc form for the graviton propagator. We begin 
by expressing the two primitive diagrams of Fig. [T]in terms of propagators 
and vertices. We next point out that the form of these expressions lends 
itself to a simple representation for x') as the sum of two tensor 

differential operators acting on structure functions. We express the BPHZ 
counterterms (the third diagram of Fig. [T]) directly in terms of their contri¬ 
butions to these structure functions. The section closes with the derivation 
of an important identity concerning the photon propagator which permits a 
great simplihcation of the more difficult, hrst diagram of Fig. [1] 

2.1 Notation and Primitive Diagrams 

Although we will cite original work, a unihed treatment can be found in 
section 5.2 of 157]. We work on the spatially hat cosmological patch of de 
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Sitter whose invariant element is, 


ds'^ = = ({^{—drf + dx ■ dx) 



( 8 ) 


where H is the Hubble constant. Note that g^y = where g^y is the 

Minkowski metric. We work in D spacetime dimensions to facilitate the 
use of dimensional regularization. Whereas the D — 1 spatial coordinates 
—oo < X* < +00 take their usual values, the conformal time g runs from 
g —)• —oo (the inhnite past) to 77 ^ 0 “ (the inhnite future). 

The vacuum polarization i[^n'^](x; x') is a bi-vector density which depends 
upon two spacetime points, x^ and x'^. In representing functions such as 
propagators which depend upon these two points, we will make extensive use 
of the de Sitter length function. 


y{x;x') 


a{g)a{g')H^ 


X —X 



(9) 


We also need the de Sitter breaking product and ratio of the two scale factors, 

u{x; x') = \n(aa') , u(x; x') = In^—j , (10) 

where a' = a{r]'). The de Sitter metric at x^ and x'^, along with products 
of derivatives of y (without the ie term) and u furnish a convenient basis for 
representing bi-tensor functions of x^ and x'^, 

di^y , d'^y , df,dly , , d'^u . ( 11 ) 


(We do not require derivatives of u(x;x') because d^v = +d^u and d'^v = 
—d'^u.) It turns out that either taking covariant derivatives of any of the hve 
derivatives m, or contracting any two of them into one another, produces 
more elements of the basis [5H1 EH]- 

The Maxwell Lagrangian is ^Max = — g and the only 

interactions we require descend from its second variation. 


S^S- 


Max 


Mp(x)Mp(x' 


= -d.d[ 




gKAgMp_ 




siD, 


X —X 


( 12 ) 


The necessary vertex functions are obtained by expanding the full metric 
around the de Sitter background as in ([2]). We can take advantage of the 
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conformal flatness of the de Sitter background to extract the 

scale factors and express the result using the notation of previous work in 
flat space background [39] and on de Sitter, with the conformally rescaled 
graviton held, in the noncovariant gauge [H], 

+ 0{k^) . (13) 


The tensor factors for the 3-point and 4-point vertices are 

yppnXap _ ^K[X^p]p y^^a)[p^K][p^X]{li ^ 


(14) 


jj-fipK.XajS'yS _ 






-|-^75^a) [p^k] [p^X] (p _|_ ^K(a^/3) [A^p] (7^5)p _|_ [A^p] (a^/3)p _|_ ^k(«^/ 3)(7^<5) [A^p]p 

-l-y^K(7^5)(a^/3)[A^p]p _|_ ^k[A^p](o^/3)(7^5)p _|_ ^K[X^p\{'y^&)(a^p)p ^ 


Parenthesized indices are symmetrized and indices enclosed in square brack¬ 
ets are anti-symmetrized, and both are normalized. 

We can express the hrst two diagrams of Fig. [Tjusing the vertices fll41ll5p . 
along with the graviton propagator i[ap/S.y] (x; x') and the photon propagator 
i[pAo-](x; x'). The leftmost diagram is formed from two 3-point vertices. 


PJTU 


{X]x') = dM i 


aP^-fS 


(x; x') 


dxd'^i 

The middle diagram contains a single 4-point vertex. 


p^a 


{x-,x') } . (16) 




4pt 


= d.d[ 




ap^-yS 


(x; x) 6 


^^x—x') 


(17) 


2.2 Representing the Tensor Structure of z[^n^](x; x') 

It can hardly escape notice that each of the two primitive diagrams fll61ll7p 
takes the form of one primed and one unprimed derivative contracted into a 
bi-tensor densityjl 


^Note that our quantity x') contains a factor of yj—g{x) x ^/—g{x') = {aa')^ 

that was not part of the symbol of the same name employed in Ref. [601 [61] . 


(x; x') = dpd'^ 


pprj^V 


(x; X 


(18) 
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which is antisymmetric on each index group and symmetric under reflections, 


pprjil/cr 

[x] x') = — 

pprjnlXT 

{x] x') = — 

pprjiau 

{x-,x') = + 

uarppp 


(19) 

(Note that because the vacuum polarization is a bi-vector density, there is 
no distinction between divergences formed with ordinary or covariant deriva¬ 
tives.) This form 01811191) is the necessary and sufficient condition for the 
vacuum polarization to be reflection symmetric and transverse, 


df,i 




= 0 = d'J 




( 20 ) 


The form 0181I19P pertains to any background geometry, not just to de 
Sitter. When the background metric and the propagators possess isometries 
great restrictions on the form of the bi-tensor x') can be imposed. 

For example, the flat space result ([5]) shows it can be reduced to the form of 
a single tensor times a scalar structure function. 


flprrii/a 
^ fiat 








( 21 ) 


This form fl^T]) is obviously more economical than acting the derivatives in 
(fT8|) and writing out all the resulting tensors. It is also more straightforward 
to employ in the quantum-corrected Maxwell equation ([1]). Not the least of 
this representation’s advantages is that the one graviton loop contribution 
to the structure function J^(Aa;^) is only quadratically divergent, as opposed 
to the quartic divergences in the primitive diagram. 

The de Sitter geometry has as many isometries as flat space so one might 
expect that a similar representation is possible in terms of just one struc¬ 
ture function. However, a second structure function is required because the 
graviton propagator breaks de Sitter invariance [33 [38] down to just the 
cosmological symmetries of homogeneity and isometry. We have chosen to 
represent the result in the same form which was first employed for the one 
loop contribution from scalar quantum electrodynamics on de Sitter back¬ 
ground [62l [63] , 


'"’T'"’] (i; x') = xF{x; x') + (i)'“'5j'”-ij'"))'''’) x G(x; x') . 

( 22 ) 

Here and henceforth, an overlined tensor represents the suppression of tem¬ 
poral components, for example, + SqSq. Our representation (l22i) 


















has a transparent physical interpretation [35] and is simpler to use even in 
cases for which a de Sitter invariant representation is possible [50] • There is 
also a straightforward procedure for passing from one representation to any 
other [6T], so nothing is lost by employing the form fl2^ . 

The primitive diagrams fll61ll7p each permit one to read off a contribution 
to the tensor x') , however, this contribution is not immediately in 

the form fl22D . The most general form consistent with homogeneity, isotropy 
and the symmetries (fT9|) can be expressed as a linear combination of the 
basis tensors formed from differentiating y{x;x') and u{x]x') [601 IGT] . 

^ x My,u,v) + D^^y DP^D'^^y x h{y,u,v) 

+D^PyDP^D'^'"yD'‘"^u x f 3 {y,u,v) + dK D^^D'^'^y D'^'^y x f 3 {y,u,v) 
+dK DP^D'^^'y D'"^u X U{y, u, v) + D^Py DP^u D'^'^yD'^'^u x U{y, u, v) , (23) 

where we dehne the f3{y,u,v) = f3{y,u, —v). (Acting on scalars as in fl23|l 
we have = gP'^{x)da, D’’^ = g'^^{x')d'p.) Only two combinations of the 
fi{y,u,v) are independent. We call these the “master structure functions” 
$(j/,M,n) and '^{y,u,v). Given the values of the various fi{y,u,v), one con¬ 
structs the master structure functions according to Table H] [5T] . 

One constructs the structure functions F{x;x) and G{x;x') of our repre¬ 
sentation according to the rules [6T|FI 


F{x;x') = {aa')^ ^xJ 

-2$ 


G{x-,x') = {aa')^~^ X 

\d- 

-l)<h -h ydy^ + 2du^ + T 


(24) 

(25) 


Here and henceforth the symbol /[/] represents the indehnite integral with 
respect to y, 

= j dy'f{y',u,v) . (26) 

Note from expressions fl241l25p that the structure function G{x; x') is less 
divergent than F{x;x'). Whereas F{x]x') must possess quadratic ultraviolet 
divergences which are determined by the flat space limit, G{x]x') is only 
logarithmically divergent, and it vanishes in the flat space limit. 

^The formulae we give contain factors of (aa')P‘~^ which were mistakenly omitted from 
equations (38-39) of [61] . 
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i 

4W-4 X (!>i{y,u,v) 

AH Ti(|/,M,n) 

1 

2{D-l)f3-2{2-y)dyh 

2(SJ - 3l)h 

2 

D{2-y)f2 + {Ay-y^)dyf2 

(2 - y)[dl - dl)h 

3 

-{D-l){f3+h) + {2-y)dy{f3+f3) 
-2dy(e'’ f3+e-^ f3) 

2dydu{e'’f3+e-'’j3)+2dyd^{-e'’f3+e-'’j3) 

-{dl-dl){h+h) 

4 

dyU 

-{D - l)dyfi + (2 - y)d^fi - 2dyduU 

5 

-Df3+2{2-y)dyh 

—4cosh(i))9y/5 

{D-l)f5-{D+l){2-y)dyf5-i‘iy-y^)dif5+2duf5 

-l2i2-y)-icosh(v)]dyduf5-4:dydv[smh{v)f3]-(d^-d^)f3 


Table 1: The contribution to master structure functions $(?/, m, v) = 

and \h(r/,M,r;) = Yl^=i from each of the coefficient functions fi{y,u,v) of 

equation (|23]). Note that f3{y,u,v) = f3{y,u, —v). 


2.3 Renormalization 

Einstein + Maxwell is not perturbatively renormalizable [Ml Eh], but we can 
still absorb the divergences using BPHZ counterterms of higher dimension in 
the standard sense of low energy effective field theory [a Emu. Our gauge 
(Ej) is covariant and we employ dimensional regularization, so it might be 
thought that only invariant counterterms are required. That turns out not 
to be true for three reasons: 

• Our interactions are time-ordered, as is apparent from the ie prescrip¬ 
tion of the de Sitter length function (Ej) which enters propagators; 

• Quantum gravitational interactions possess two derivatives, which al¬ 
lows the noninvariant ordering to contaminate expression fllbp : and 

• The coincidence limit of the graviton propagator diverges in de Sitter 
background instead of vanishing like it does in flat space [5U] . 

We are loath to change the ordering prescription because it is so embedded 
in the Schwinger-Keldysh formalism [Ml EZl EHl EHl EOl [Ml UM EHl [M] that 
must be employed to make the nonlocal part of equation ([T]) both real and 
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causal. The only alternative is to permit the same noninvariant counterterm 
that was necessary when using the noncovariant gauge j41j . 


+CH^F,,Fp,g^Pg’^^^ + ACH^F,,FMg^^g^^V^ • (27) 


The term proportional to C is the de Sitter specialization of three countert¬ 
erms for which the two field strengths are contracted into g^^g^^R, 
and . The term proportional to AC is the noninvariant counterterm 
by virtue of its indices running only over space. After some tedious tensor 
algebra one finds that AC contributes to F{x;x') and G{x]x') as [H] , 



AF = 4a^-^ !,\c-i3D-8)C^H^-—d^ + ^ 

{i J aa' 

l\G = ia°-‘\l^C - (D-6 )c}h\5°(x-x') . 


(29) 


Although the AC counterterm is strictly only required to renormalize diver¬ 
gences from the time-ordering we will also employ it to simplify G(x; x'). The 
coefficient C^ is fixed by the fiat space limit [39] in terms of the quantities 
C 2 and C'o(&) defined in expressions ([MZD) 




C 2 + Co{h) . (30) 


647rf {D-l){D-2)\D-3){D-A) 


Note that the spin 2 contribution to C 4 , is finite by virtue of the factor of 
D — 4 in expression (jHI) . 

Those accustomed to modern techniques of renormalization in covariant 
gauges sometimes hnd the appearance of noninvariant counterterms to be 
disconcerting. However, it is important to realize that they pose no problem 
of principle. The divergent part of the counterterm is of course hxed by the 
primitive divergences it is to remove, and the hnite part can be determined 
to enforce physical symmetries. (In our case the focus on late times obviates 
the need for this as long as the finite part of the counterterm is assumed 
to be of order one.) The procedure is explained in older standard texts on 
quantum field theory, for example m- And it is important to recognize that 
many of the classic computations of quantum electrodynamics were in fact 
performed using noncovariant gauges [7B] . 
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2.4 Reducing the Photon Propagator Term 

Our photon propagator is defined in exact Lorentz gauge, 


DO 




= 0 = D' 




(x; x') 


( 31 ) 


An early result m for this propagator contains a small error which was 
corrected 20 years later |^0 The most useful form for our purposes was 
obtained in a study of the graviton propagator [37] . 


oA„ 


x;x) = 


2 m L 




In'-D'^D' 


dad'py X iABBB{x;x') 


(32) 

Here the BBB-tjpe propagator is defined by inverting the 3rd power of the 
scalar d’Alembertian with mass M| = {D — 2)H^, 


□ - {D-2)m iABBBix;x') = iAbb^x-.x') , 
□ — {D — 2)m iABB{x]x') = iAB{x;x') , 

□ - {D-2)m]iAB{x;x') = 


i5^{x — x') 


V-s(x) 


(33) 

(34) 

(36) 


Because M| is strictly positive, iAbbB) ^Abb and iAb are all de Sitter 
invariant functions of y{x-,x') [52] whose precise form can be found in [57] . 
We will give the expansion for fA^ = B{y) in Appendix A. We also give 
there the expansion for the propagator iAc{x;x') = C{y) of a scalar with 
mass M| = 2{D — 3)H‘^. An important relation exists between them which 
we will use many times [83], 


2C'iy) = {2-y)B\y) - Biy) . 


(36) 


The differential operator used to construct the photon propagator flHTD 
has the key property of transversality [57], 


DP 


6 ^n-D^Dr 


= 0 = 






(37) 


^There have been some recent false claims about this in the mathematical physics 
literature m so it is important to note that one really does need to employ the corrected 
propagator in computing standard things such as the effective potential of scalar quantum 
electrodynamics [SOllHl]. Using the uncorrected propagator would not even recover the 
famous Coleman-Weinberg potential [82] in the flat space limit. 
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Exploiting this, the de Sitter invariance of iAsBB and relation 
straightforward to show 


6 ^n-D^Dp 


□' - dc^d'py x iAbbb (x; x') 

[dp&py X iAbb{x-x') 
The next step is to act the remaining primed transverse projector, 




In'-D'^D' 


do^d'py X iA 


BBB 


= dpd'^y X iAb - H'^dpd'A i 2 -y)iABB - (D-S)! 


iA 


BB 


it is 


(38) 


(39) 


The hnal step is based on the fact that the 3-point vertex factors in 
expression flTB]) inherit an anti-symmetry from the Maxwell held strength, 

_ _YU(j>ea'y5 


_ _ -^p,XKpaP 


This anti-symmetry can obviously be communicated to the differentiated 
photon propagator in flTBl) . 


DxD't 




{x; x') 




p] Ao-]] 


[x] x') . 


(41) 


Here the single square brackets indicate anti-symmetrization under A O p 
while the double square brackets indicate anti-symmetrization under (j) ■H- a. 
Because the double covariant derivative of a scalar is symmetric we can use 
(l39|) to conclude. 






(x; x') 


1 


9[pd'[[^yix;x') X dx]d'^]]B 



(42) 


Similar reductions of the original photon propagator structure functions down 
to the iAB propagator have also been noted in explicit two loop computations 
involving the very different interactions of scalar quantum electrodynamics 


Identity fl42]) allows us to re-express the 3-point diagram flT^ as. 


P.JTU 


K 


2 H^ 




a/sA-fS 


dpd'^y X dxd'.B 


(43) 

This is as far as we can get without exploiting the explicit form of the graviton 
propagator. 
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3 Spin 2 Contributions 

The purpose of this section is to work out the contributions to the renormal¬ 
ized structure functions from the spin two part of the graviton propagator. 
We begin by describing how this part of the propagator is expressed. We next 
work out the contributions from the 4-point diagram and the local part of the 
3-point diagram which comes from the delta function part of the doubly dif¬ 
ferentiated photon propagator. A much more involved analysis is necessary 
to work out the nonlocal contributions to the 3-point diagram from the de 
Sitter breaking part of the propagator and from the de Sitter invariant part. 
The section closes by giving the combined results for F{x;x') and G{x;x'). 


3.1 Spin 2 Part of the Graviton Propagator 

The spin two part of the graviton propagator is transverse and traceless and 
takes a form analogous to the transverse photon propagator fl32|) [371138] . 




2 /D-2 




H^\D-3 


X p..,?^(x) X py 




pa 


X 


X 


dad'y{x-,x') X di3&sy{x-,x') x i/^AAABB{x]x') . (44) 


Here P is a 4th order differential operator which is transverse and 
traceless on both index groups |37l 138] . 


p 1/^-3 

2\D-2 




U-DH^ 


D-2 

D-l 

D{p,Dy) 

D-l 


n+2{D-i)y 


n-2y 


+ 


u+y 




VD-l 


yn+2y 


gO-fi _ 9p 


D-l 


n+2{D-l)y]D^^Df^^ (45) 


The spin two structure function iA^AABB^x] x') is constructed by inverting 
the kinetic operator for a massless scalar three times and inverting the kinetic 
operator for an M| = {D — 2)H^ scalar twice. The order in which these 
inversions are performed is irrelevant but we find it convenient to alternate, 


niAAAABBix;x') 



{D-2)y 


iA 


AABB 


[X]X') 


iAAABB{x]x') , 
iAAAB{x-,x') , 


(46) 

(47) 
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niAAAB{x]x') 

= iAAB{x]x') , 

(48) 

□ - {D-2)H‘^ 

iAAB{x;x') 

= iAA{x]x') , 

(49) 


□tA/i(x; x') 

i5^{x—x') 

■ 

(50) 


Scalars with M| < 0 inevitably break de Sitter invariance |86l |87l EH EH 
[59] so each of the integrated propagators in relations (14611501) does as welllfl 
That is, we can express each of the integrated propagators as the sum of a 
de Sitter invariant function of y{x;x') plus a de Sitter breaking term, 

iAAAABsix; x') = x') + i^AAABBi^', • (51) 


Explicit forms for both terms can be found in |95]. We refer to the “de Sitter 
breaking” and “de Sitter invariant” part of the spin two graviton propagator 
as derived from and respectively. The de Sitter breaking 

part is [95] . 


A 2brk 


{x]x') = {aa'Y 


D-1 


VfiuVpa 


X k 


ln(4aa')+y42 , (52) 


where we recall that and the constants k and A 2 are. 


k = 


HD -2 r{D-l) 
(4^)T r(f) 


(53) 


(Here Piz) = ^ln[r(z)] is the digamma function.) The de Sitter invariant 
part is much more complicated. The graviton analog of relation fl38|) is [95] - 


A 2inv 
pi'^pa 


[x] x') = —H 


-4 


D-2 

D-3 




X 


d^d’yxd^d'gy y n-{D-2)H‘^ 


. (64) 


®The mathematical physics literature contains claims that there is no need for a de 
Sitter breaking part isniiiiiiiij- Morrison |93j has shown that constructions which purport 
to give a de Sitter invariant propagator differ from ours in two ways: (1) the propagator 
for a scalar with general mass-squared Mg must be considered as both de Sitter invariant 
and well defined for all M|, except for simple poles at M| = —N{N + D — 1)H^ with N = 
0,1,2,..., and (2) it must be accepted that, for constructing the graviton propagator, an 
arbitrary constant can be added to equation (ISOl) . Both of these deviations are illegitimate, 
resulting in formal solutions to the propagator equation which are not true propagators 
in the sense of being the expectation values, in the presence of positive-normed states, of 
the time-ordered product of two graviton field operators [53] . 
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i 

My) 

1 

Iwlw - (D-i)\J(y) 

2 

I[-{2-y){^jy + (D-l)gJ + {D-2){2-y)J' - {D-l){D-2)J] 

3 

/21-i§J + 2(g5f)J" - \(D-2)(2-y)J' + \(D-2)(D-l)J] 

4 

ThwXy) + r-[-\(2-y)(^J)' + |(I>- 2 )§ J] 

5 

D-i^M{y) + \{‘^ y)‘^^2J{y) \{D 3)(2 y)l[^2J] 


Table 2: The five functions Ji{y) in relation fl56|) for the de Sitter invariant 
part of the spin two graviton propagator. The function J{y) is given in 
expression (155|) and /[/] stands for the indefinite integral of f{y) as in fl26l) . 


In Appendix B we show that acting the derivatives gives, 

— — a;') = x') + constant = J{y) . (55) 

The analog of relation fl5^ is, 


A 2inv 


+H'‘ g;,„D'^D'^+g'^Di,D^ Ja(g) + H‘‘g^^g'^ J^{g) , ( 56 ) 


where Table |2] gives the functions Ji{y). 

Relation fl56|) is remarkably similar to relation fl3^ . In particular, the 
tensor structure of the first term is provided by second derivatives of y{x-, x'), 
and the function Ji{y) is 2 A™''/ 2 if^ plus a constant. The remaining terms 
are all gradients and/or traces. Unfortunately for us, neither gradients nor 
traces drop out of the vacuum polarization the way the analogous photon 
gradient terms did in expression (H2|) . This leaves no alternative but to act 
the various covariant derivatives in expression fl56p and express the result as 
a linear combination of the de Sitter invariant bi-tensors given in Table [21 


A 2inv 


E 

2 = 1 


pu I pa 


{x-,x') X Ki{y) 


(57) 
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i 

[f.uT;j{x-, x') 

1 

dpd[py{x; x') d'^^d^y{x-, x') 

2 

d(py{x; x') d^)d[py{x; x') d'^^yi^x] x') 

3 

dpy{x] x') dyy{x] x') 9(,|/(a;; x') d'^y{x] x') 

4 

H‘^[ 9 ^lu{.x) d'^yix] x') d'^y{x-, x') + d^y^x] x') d„y{x; x') gpa{.x')] 

5 

H^9tiu{x)gp„{x') 


Table 3: The de Sitter invariant basis bi-tensors in relation (jST]). As always, 
indices enclosed in parentheses are symmetrized. 

The coefficient functions Ki{y) are given in Table 01 It is worth noting that 
tracelessness implies two relations among the Ki{y), 

4Ki + {4y-y^)K^ + DK, = 0, (58) 

-K, + {2-y)K2 + iAy-y^)Ks + DK, = 0. (59) 


3.2 The 4-Point Diagram 

The simplest diagram is the middle one of Fig. [1] The spin two contribution 
to it comes from substituting the spin two part of the graviton propagator 
into expression (HU. Because the spin two part of the graviton propagator 
is traceless at each point, we can drop terms in the 4-point vertex 
which contain either or Many of the other terms are also related 
when contracted into i[aj 3 ^.y&\{x]x) so that there are only three distinct con¬ 
tributions. 




4pt2 


X 


{x-,x') = dpd'^l —K^a^ H6^{x—x') 




q/3 A 2 
^ 0/3 


+ 2aH 


.(60) 
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i 

Ki{y) 

K^{y) 

1 

Ji + J'2 + 2J'^ 

+ i{D-2){Cj)+4(g52)J'' 

2 

JP+AJ'P 

+8(§5f)^"'-(^?-2){2-3/)J"+(D-2)DJ' 

3 

jnn 

-5(;PT"+2(+f)J"" 

-i(D-2)(2-j/)J'" + i(D-2)(D+l)J'' 

4 

-j^+|(2-s)4]"+j: 

flir (ft -^)"+2 (gef) (2 - j,) J'" - J" 

5 

~ ( 2 -i;) J 2 + ( 2 -i/) [ ( 2 -^) J 3 ] ' 
+ 2 ( 2 —y) J^+J 5 

--+T(^T+-+T( 2 -y){-^T' 

- 2 (+f)(ST+ 8 (+f)-^" 


Table 4: The coefficient functions Ki{y) expressed first as derivatives of the 
functions Ji{y) from Table [2l then in terms of the function J{y) dehned in 
expression fl55ll . 


Recall that indices enclosed in square brackets are anti-symmetrized, and 
that the double square brackets in the hnal term of fl60l) serves to distinguish 
the anti-symmetrization on i/ -f-)- a from that y p. 

The coincidence limit of the graviton propagator takes the form [95], 




s%s%+sv\-- 




+ 




X iAi 

2 —a0— 


X i/\2{x) . (61) 


Hence the three terms on the second line of 


are, 




ay A 2 
^a 0 


{D-l){D-2)iAi 

+ {D-2){D + l)iA 2 {x)\ X pp^’^p^^P , (62) 


2 a^i 


^[u ^2a-]p 


= -2 


D + 2\ . 

“d” 


iAi X 


D + l\ . 
D-l 


iA2{x) X , (63) 




4(D-l)(Zl + 2).^ , 2(Z1-2)(D + 1).^ , , 

—H--—;- iA2(x) 


D 


D-l 
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X . (64) 

In deriving the last of these relations we have used, 

• (65) 

Expressions 0601) and 062tl64p imply that the 4-point contribution already 
takes the form described in subsection 2.2, 

= dpd'^!^27^^^^''r]^^PxF2a{x-x') + 2fJ^^^^^ . ( 66 ) 

The two structure functions are. 


F 2 a{x;x') = {x — x') 


{D + 2){D-l){D-A) . 


{D+2){D-2) . 


D 


iAi 


iAi 
AD 

{D-h){D-2){D + l) . 


4(B-1) 


/O / ^\ 2 D—A • / /\ 

G 2 a{X',X) = Ka to (x — x) 


(D-3)(D + 1)^. 


D-1 


iA 2 (x) 


iA2(x) [> , (67) 

( 68 ) 


From the previous subsection we hnd the constant iAi and the time depen¬ 
dent function iA 2 {x) to be, 

4{D-2)D{D+1) 


iAi = 


iAo = k 


{D-1) 

2 ln(2a) -|- A 2 


i7V"(0) , 


= k 


2 ln(a) -|- A 2 


(69) 

(70) 


Comparison with expressions 0281I29P suggests that we choose the “2a” con¬ 
tributions to the C and AC counterterms to be, 

K^iAi {D + 2){D^-9D + 12) iFkA2 {D-3){D-2){D+ 1) 


Co„ = 


772 


1677 


772 


16(77-1) 


AC2a = 


n^kA2 (77 - 3)(77 + l) 


772 4(77-1) 

The renormalized 4-point contributions to the structure functions are, 

5^2772 


FTaiDx') = 

G'^27{x-,x') = 


247r2 

5^2772 


127r2 


ln{a)i6'^{x—x') , 
\\i{a)i5^{x—x') . 


, (71) 

(72) 

(73) 

(74) 
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3.3 General Form of the 3-3 Diagram 

The 3-3 diagram is the leftmost part of Fig. [T]and is by far the most difficult 
to evaluate. The hrst step is to substitute the 3-point vertex factor ffTTj) into 
expression fH3|) . Because the spin two part of the graviton propagator is 
traceless we retain only the second terms, 

ylipnXa^ _^ ^k][p^ yuaetfj'yS _^ 


It is desirable to expand out the anti-symmetrizations over p A and a (j), 
and also to re-label the external derivatives from d^^dg to dpd'^, 


PJTU 

■'^J-3pt2 


(x; x') = dp&„ 


2 K^{aa 

IP 


/\D r 


[p 

/3 


{x; x') 


DP^D'"'^\ X D^D'^B 


X DP^D'^B-PP^D'^y x + x PP^D'^'^^B 


(76) 


In expression fl76|l we have employed single and double square brackets to 
distinguish the anti-symmetrizations over p -f-)- p (single) from that over 
u -H- a (double). 

Recall from subsection 2.2 that the structure functions can be computed 
directly from the portion of fl76|) within the curly brackets. It is useful to 
give each of the four terms its own symbol, 


pprp 


ua 

2 


pprpl/cr 

2-^2 


Pprpua 

3-‘-2 


[X] X ) 

{x; x') 
{x] x') 


fip/jiucr 

4-^2 


{x] x) 


2 p(aa')^ 

IP 

2 p(aa')^ 

IP 


X t 


[p ^2[[v 



2 p{aa')^ 

IP 



2 p{aa')^ 

IP 


X i 



X X P^P'^B , (77) 

X P^P'^'^^y X PP^P'^B , (78) 
X pP^P'^y X , (79) 

X P^P'% X . (80) 


Our notation is that the left hand subscript denotes which of the four permu¬ 
tations is intended, while the right hand subscript “2” indicates that these 
are all contributions from the spin two part of the graviton propagator. 

The next step is to act the derivatives on i?(p) [96] . 


P^P'^B{y) = 


-^0 "0 


a 


D+2 


■ + D^yP'Py 
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It is natural to distinguish the “2b” local delta function terms (subsection 
3.4) from the nonlocal terms. We further distinguish the (nonlocal terms 
between the “2c” ones from the de Sitter breaking part of the graviton prop¬ 
agator (subsection 3.5) and the “2d” ones from the de Sitter invariant part 
(subsection 3.6). 


3.4 Local Contributions from the 3-3 Diagram 

The “2b” contributions come from the replacement in expressions fl771l80p FI 

2 H^ 


DPD""y X D^D'^B 


,0+4 




(82) 


The delta function re-introduces the coincident graviton propagator expres¬ 
sion fIM]) . With replacement (15^ the four contractions in fl76D give, 


D-2 

D 




jJ^prjii/a 
2-^ 2b 


flprpi/a 

3-^26 


p,prpi/a 

4-^26 






(83) 

(84) 

(85) 


= An'^a^ '^i6^{x — x') 




D D-1 

We can read off the structure functions using the relation, 


( 86 ) 


(87) 


®Note that the analytic continuation of Barvinsky and Vilkovisky would give a 
covariant result instead, 


DPO^'y X D°‘D'^B 


pa a 0 iS^jx-x') 
D ^ ^ ^ 


Because physics is ultimately based on the Minkowski signature of our expression (jH]) we 
feel it is safer to work with (I82|) . even though it entails a noncovariant counterterm. Some 
of the problems which can arise from analytic continuation are explained in |94) . 
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The results are, 


F 2 b = n^a^ [x — x') 

G2b = K^a^~'^i6^{x — x') 


{D+4){D-2) , {D + l){D-2) , 

-- Z/A 2 [X) 


D " {D-l) 

(D+l){D-2)^ 


D 


■ iAi 


{D-l) 


D^ 2 {x) 


( 88 ) 

(89) 


The constant iAi is divergent, which is why this computation requires 
a noninvariant counterterm. Comparing expressions fl881l89p with (12811291) 
suggests that we take the “2b” contributions to C and AC as. 


- K^i^i{D+A){D-2) K'^kA2{D + l){D-2) 

~ Id w Ip 4 ( 0 - 1 ) ’ 

^ kHAi{D^-8 ) PkA2{D + l){D-2) 

~ AD ^ m 4(0-1) 

Our hnal results for the renormalized structure functions are, 

= —^^\n{a)i6'^{x-x'), 

G3“(i;i') =- \B(a)ii‘(x-x') . 


(90) 

(91) 


(92) 

(93) 


Note that the noncovariant divergence in expression (1891) is exactly cancelled 
by the noncovariant counterterm (I?I1) . so there are no spurious finite terms. 


3.5 Nonlocal de Sitter Breaking 3-3 Contributions 


The “2c” contributions derive from making the following replacements in 
expressions (I771l8nj) . 


i 


/\2v 

13 


{x; x') 


D^D'^B{y) 




2 


X k 


ln(4aa') + A2 


D^D'^y X B'{y) + D^y D'^y x B"{y) . 


. (94) 
(95) 


Recall that an overlined tensor indicates the suppression of its temporal com¬ 
ponents, — ^o^a- These overlined tensors in expression (l9T)) can all 

be represented using the standard basis described in section 2.2 [95l [98] , 


= 


aa 


-Df^D'y+Df^y D'u + D^^u D'y + {2-y)D>^u D'u 


(96) 


22 




















fi 

nT,7] 

r/ip 

[2+3-^ 2c J 


[ppy..] 

fl 

(D-2)(D+1) 

2(D-1) 

yD-i) 

{D-2)(D+l) 

2{D-1) 

(D-3)(D+1) 

D-l 

f2 

0 

0 

0 

0 

/s 

(D-2)(D+1) 

2(D-1) 

yD-i) 

(D-2)(D+l) 

2{D-1) 

(D-2.){D+1) 

D-l 

/s 

(D-2)(D+1) 

2{D-1) 

\D-l) 

{D-2){D+l) 

2{D-1) 

{D-2.){D+1) 

D-l 

u 

'h'w'-ir'p V) 

‘i{‘2-y) 

+ 0 ^ cosh(t)) 

‘°2TT°r’(2 V) 

y) 

+ 0 ^ cosh(D) 

fb 

0 

/D-3N 

^D-l) 

0 

/D-3N 

^D-l) 


Table 5: The contribution proportional to B'{y) from each permutation type 
to the coefficient functions fi{y,u,v) which were defined in expression (|23|) . 
Each contribution should be multiplied by B'{y) x + 21n(2) + A 2 ]. 


The hrst step is to substitute relations fl9T|l . fl95|) and fl^ into the stan¬ 
dard permutations fl771l80p and read off the contributions to the hve coef- 
hcient functions fi{y,u,v) which were dehned in expression ([23]). Relation 
fl^ contains a term proportional to B'{y), whose contributions to each /j is 
given in Table O The B"{y) contributions are listed in Table [HI 

The next step is to compute the master structure functions u, v) and 
'^{y,u,v) according to the rules which were originally derived in Ref. [61] 
and which are summarized in Table [Hof section 2.2. The intermediate results 
can be substatially simplihed using the R-type propagator equation, 

{Ay-y^)B"{y) + D{2-y)B\y) - {D-2)B{y) = 0 . (97) 


The final result for the “2c” contribution to <h(|/,M,n) is, 


<h 2 c(|/, M, v) = K^H^k ( ^ -1) u + 2 ln{2) + A 2 


xM 


{2-y)B cosh(n)-8R"'-(R-3) (2-|/)R M98) 


The “2c” contribution to T(?/,M,n) is more complicated because it involves 
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fi 

[TUTl 

r/ip rppfjl 

[2+3-^ 2c J 


r^2T] 

fl 

(2-y) 

—2 cosh(4;) 

0 

0 

i2-y) 

—2 cosh(i;) 

f2 

1.D-3N 

2\D-1^ 

\D-l) 

(D-2)[D+1) 

2(0-1) 

(D2-2D-7) 

2(0-1) 

h 

\{^){2-y) 

-2(g5f)e’'-e-“ 

|(S5f)(2-y) 

2 (s—'v 

D-1^ ^ 

0 

2(%E^)(2-y) 
-2(g5|)e^-2e—" 

fs 

U^X2-y) 

-e'^-2(§^e-'" 

^(S5f)(2-y) 

-e- + T^,e- 

0 

2(S5T)(2-y) 

-2e”-2(S5|)e-'" 

u 

2(S-i)+itn-i)(2 y? 

-4(§5f)(2-?/) cosh(i;) 

0 

0 

2(o-l)+2(D-l)i2 yd 
-4(§5f)(2-y) cosh(i;) 

h 

0 

-(^)(2-y) 

-4(-§5f)cosh(D) 


2(0-1) y> 

+4(§5f)cosh(i;) 


Table 6: The contribution proportional to B"{y) from each permutation type 
to the coefficient functions fi{y,u,v) which were dehned in expression fl2^ . 
Each contribution should be multiplied by B"{y) x + 21n(2) + A 2 ]. 


derivatives with respect to u, 


'^ 2 c{.y,u,v) = 


D-l 

+ (D-3)(D-2)(D-1) (2-2/)i? 


2{D^ — ?>D+A) —{2 — y)BcoAi{v) + 2B 


+ 2K^H‘^k 


D-2 

D-l 


u + 2 ln(2) + ^2 


X 


8 B’+ {D-l){2-y)B cosh(u) - 4 (2-j/)5" 1.(99) 


The structure function F{x;x') is constructed by integrating ^{y,u,v) 
with respect to y according to the expression This is simple to do 

because each term in fl98p is either a 2nd or 3rd derivative with respect to y, 


F2c{x]x') = {D — 2)K^H‘^k{aa')^ ^ M + 21n(2) + 742 


X 


{2-y)B cosh{v) + 8B" + {D-3) {2-y)B WlOO) 
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The other structure function follows from substituting 09811991) in fl2SD . 


G2c{,X]x') = 

'D-2 


{2 — y)B cosh(u) — 2B' 


D-1 


2{D ^)D 2^2 
D-l ^ ’ 

K^H‘^k{aaY~‘^\u + 2\n{2) + A 2 \ IQB - 4:{D-l)yB cosh(u) 


165 + 4(5-3)2/5 - 4(5-3)25' + (5-3)(5-l)(j/5)' ^.(101) 


These results 01001I101|1 are valid in D dimensions. Because the quantum- 
corrected Maxwell equation ([1]) involves integrals of F{x] x') and G(a;; x') 
with respect to we can set 5 = 4 for any part of the structure functions 
which diverges less strongly than l/{x — x')^ as x' —)■ x. Expression O204p of 
Appendix A gives the expansion of B{y), 

uD-2 _ 1 / s 

B{y)= n J +o({D-A)y^) . (102) 

Because y = aa'H^Ax"^, it might seem that the terms in fllOOtllOl]) which 
involve B"{y) harbor quadratic divergences, while those involving yB"{y) 
and B'{y) diverge logarithmically. In fact all divergences cancel. To see this 
first note that the factor of cosh(u) can be rewritten as, 


cosh(u) = 


a a 
+ - 


= 1-1- ^aa'H^{ri — ri'y = 1 -|- ^aa'H'^Arf . (103) 


2 La' a J 

This demonstrates that the quadratic divergences of F2c and G 2 c cancel, 

(2-2/)5j''cosh(u) + 85" + {D-3)\{2-y)B]' = -Aaa'B” 

+4(?/5)" + 2(5-3)5' + 2aa'H^Ay\yB)" - {D-3){yBy ,(104) 
165 — 4(5 —l)j/5 cosh(u) — 165-|-4(5 —3)j/5 —4(5 — 3)25' 

+ (5-3)(5-l)(|/5)' = 3aa'H^Ay^B" - 3{D-2){yB)" - 4(5-3)25' 

-2{D-l)aa'H‘^Aif{yB)'' + (5-3)(5-l)(?/5)' .(105) 


Cancelling the logarithmic divergences is more subtle. Relations fll02p 
and (I103p suffice for the u-independent part of G 2 c, 


{2 — y)B cosh(u) — 25' 


aa'H'^ATfB' — (?/5)'cosh(u) 


1 Aif 


aa'Ax'^ 

(106) 
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In the M-dependent parts we can ignore {yB)' and Ar]‘^{yB)" which are both 
integrable and vanish in D = A dimensions. Further, we need only retain the 
hrst term of fll02p in evaluating the combinations of Arj'^B", {yB)'' and B' 
which appear in expressions (11041) and (11051) . 


-Aaa'H^ATfB" + A{yB)" + 2{D-3)B' 


H^-^T{A) f 2Daa'H^Ay^ 2 

y2 


An 2 


-+i 

y^ 


Saa'H'^Ay^B" - 8{D-2){yB)" - A{D-3)^B' 


TD-2-rfD^ 

~A 

An 2 


H^-^T{A) f ADaa'H^Arj^ ^ A 

{ + — + ••• 


y- 


y- 


,(107) 


. (108) 


Expressions (I107p and (llOSp are proportional to the same function which can 
be reduced to a form that is integrable in D = 4, 


Daa'H‘^Ari^ 1 _ If DArj^ 1 1 

l/T+i yn {H‘^aa')^\^^^^‘^ Ax^ j 

1 7 ^ An^i6^{x — x')^ 

~ {D-2){H^aa')^{ ^ r(f-l) J 


(109) 

.(110) 


After some simplihcations our hnal unregulated forms are. 


F-^{x;x') 

GlT{x-,x') 




In 


lOvr^ - 
( H^aa' 


/ aa' \ 1 1 

(-) + 3-H 




( 111 ) 


247r4 


{d^ + V^)\n{H^Ax^) 




3.6 Nonlocal de Sitter Invariant 3-3 Contributions 


The “2d” contributions derive from making the following replacements in 
expressions (I771I80I) . 


fj. \2v 

OL^ y 


(x;x') J2h'T‘Ax;x')xKfy), 


1=1 


D'^D'^B{y) 


D'^D'^ X B'{y) + D^y D'^y x B"{y) . 




(113) 

(114) 
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i 

(A0i)i(|/) 

{A(j)2)i{y) 

1 

-8{D-1)B' + 2{D-2){2-y)C' 

ADB" 

2 

iD-2)iAy-y^)C' 

A{D-2)B' - 2{D-2){2-y)C' 

3 

0 

-A{D-2){Ay-y^)C' 

4 

0 

-IQ{D-2)C’ 

5 

-8B' 

8B" 


Table 7: The functions (A0i)j(|/) and (A02)i(2/) defined in equation (11151) . 


Recall that the basis tensors [^yTp„]{x] x') are listed in Tableland the func¬ 
tions Ki{y) are given in Table HI Because everything is de Sitter invariant, 
the only coefficient functions /* from expression fl2^ which occur are fi{y) 
and f 2 {y)- It is simplest to extract a factor of and to report results 

from each of the basis tensors x') in the form, 

5 5 

/,(!/) = xA'ite) , 

i=l i=l 

(115) 

Our results for the functions and (A02)i(?/) are given in Table [3 

Relations fl36|) and fl97|) are helpful in simplifying these expressions. 

The next step is to construct the master structure functions. Only ^{y) 
is nonzero because everything is de Sitter invariant. We have found it use¬ 
ful to extract a factor of and to distinguish the terms which are 

proportional to Ki{y) from those which are proportional to K-{y), 


2=1 


^2diy) = UiD-l){A<P,),-2{2-y){A<P,)[ + D{2-y){A<P,), 


+ (%-2/^)(A02)' 


X Ki -\- 


-2(2-|/)(A0i), + (%-y2)(A02). xKi .(116) 


Our results for the coefficients of each Ki(y) and Kl{y) are reported in Table 
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i 

2{D-l){A0i)i-2{2-i/)(A</.i)' 
+D(2-y)(A</)2)i + (4i/-i/2)(A</)2)' 

-2(2-2/)(A0i)i + (4|/-|/2)(A02)i 

1 

-8{D-l)B' + 8{2-y)C' 

8{2-y)B'-8{D + 2)C'+A{Ay-y‘^)C' 

2 

4:(D+2)(2-y)B' 

-32C"+12(4y-y2)C" 

A{Ay-y^)B'-A{2-y){Ay-y‘^)C' 

3 

lQ{Ay-y‘^)B'-lQ{2-y){Ay-y^)C' 

-A{Ay-y^fC’ 

4 

UB' -82{2-y)C' 

-16(4|/-|/2)C" 

5 

0 

-16C" 


Table 8: The coefficients of Ki{y) and Kl{y) in equation (11161) . An overall 
factor of {D — 2) has been extracted from each term. 


[HI Substituting the Ki{y) from Table 01 and adding everything up gives, 


xlSB' 


D- 
\H^ 


j) +{D-2){2-y)J"'-{D-2){D-l)J'‘ 


(n 


/□2 


+ (D-2)C'[i^—j) -2{2-y)[^—j) +D[—^J) +S{D-2)J"'\ |.(117) 

The expansions of B\y) and C'{y) are given in expressions (1206112071) 
of Appendix A. The function J{y) was defined in expression fl55l) and its 
expansion is given in expression fl229p of Appendix B. When the various 
factors in equation flll7|) are combined the result is, 

^2^2D-2r2('D^ (T)-2)(T) + 1) r {D-A){D + 2) 


^2d{.y) = 




t 


(C-l) 

(3T)2 26T) + 52)T)2 


Ay 


D 


■ln(!) 


+ 


(118) 


A8{D-2)y^-^ y‘^{A-yy y‘^iA-y) 

For this case the master structure function A/{y,u,v) vanishes so we have. 


F 2 dix; x') 
G2d(x;x') 


-2(aa 


r\D-2 : 


^2d{y) 


{aa') 


l\D-2 t2 


(r>-l)'S> 2 j(!/) + !/42j(!/) 


(119) 

( 120 ) 


































Note that the term in flllSp drops out of expression fll20p . 

Substituting flllSp into (11191) . performing the integration and recalling 
that y = aa'H^Ax^ gives, 


F2dix;x') 


K2r2(f) {D-4){D-2){D + l){D + 2) 
327r^ {D-iyaa'Ax^^-^ 


{3D^-26D + 52)D^{D + 1) 
3847r^ {D-2){D-l)Ax^^-^ 

5K^H%aa'f f 1 , [1 1 1 

+ 144^- 




.( 121 ) 


where we dehne the function C{y) as, 


C{y) = L^d) +ln(l 



( 122 ) 


The dilogarithm function Li 2 (^) is dehned in (12281) . 

Renormalization is accomplished by hrst localizing the ultraviolet diver¬ 
gence by a combination of partial integration, adding zero in the form of a 
delta function identity, and taking H = 4 in the hnite, integrable remainder 
[96], 




2{D-3){D-4) lAx^^-^i ’ 


+ 




H6^{x — x') 

2(T)-3)(T)-4)r(f-l) ' 2(T)-3)(T)-4) LAa;2^-6 


Att^ S6^{x — x') 

2{D-3){D-A)r{A-l) 


d"^ rln(p^Ax^) 
4 - Ax"^ 




D-4 T 


Ax^~‘^- 


(123) 

,(124) 

(125) 


We also employ the relation. 


1 Att^‘^ d^i6^{x—x') 

Aa;2^-2 ^ 4(T)-2)2(T)-3)(T)-4)r(f-l) 


rln(p,^Aa;^)‘ 
A^2 . 


(126) 


The localized ultraviolet divergence is absorbed by a counterterm of the form 
fl28P with the and C coefficients. 


KV^-^r(f) {D + l){D + 2) 
647r# {D-mD-3) ’ 


(127) 
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967r? 


D‘^{D + l) 


(D-l)(D-3)(D-4) 

{?>D-lA)D^{D + l) 3(3D-8)(D + l)(D + 2) 


+ 


4(D-l)(D-3) 
The final renormalized result is, 

c:^2rT2 

FlT{^-,x') = Ha)i5\x-x') + 


2(T)-1)2(T)-3) 
ln(/i2Ax^)' 


.(128) 


327r2 


bK^H\aa') 




Aa;2 


4£fe) 


rl 

y 


1447r^ 

The second structure function is, 

,,(D-i)(D + l)(D+i) 




G2d{x-, x) = 




6471--D 


{D-iyAx^^-^ 


bn'^H^aa') 

72^^4 


-{l-y)C{y) + 


A-y 




,(130) 


We go through the same procedure of localization as before, and add a coun¬ 
terterm of the form (|2^ to produce the final renormalized result, 

,/\2 


Gl7{x;x') = 


bt7H^{aa'f f 1 

7^ 


-(l-y)£(y) + 


-y 


In 


(131) 


3.7 The Full Spin 2 Contribution 


It remains to combine the various spin two contributions to the two structure 
functions. Our results for F 2 {x-,x') are relations fl73|) . fl9^ . fillip and fll29p . 


oc:^2 tt 2 

F 2 {x]x') = ——— \n.{a)i5'^{x—x' 


k^H 


2 r 


In 


/aa'\ 1 ^ ■ 


727r2 Ibvr'^ 

5k^H^ ^w\n{y^Ax^)\ 5K^H%aa'y (C{y) 2(2-|/)ln(f) 


Aa;2 


1447r4 




{■ 


Aa;2 


/ 1447r"^ \ 


+ 


4J/-J/2 


(132) 


Our results for the other structure function are given in equations fl74p . fl9^ . 
(mD and ([m]). 


G2{x;x') = 
K^H'^aa' 


47r2 


\n{a)i6'^{x — x') -|- 


k‘^H 

24^ 


2 r 


aa 


In — ]+--2'j 


Ax2 


967r4 


{d^+V^)\n{H^Ax^- 


5K^H%aa'f jil-y)Ciy) , (|/-3)ln(|) 


727r4 


-y 


( 133 ) 


Note that there can also be arbitrary finite contributions of the form 
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4 Spin 0 Contributions 

This section follows the analysis of Section |3] applied to spin zero part of 
the graviton propagator. Its purpose is to work out the contributions to the 
renormalized vacuum polarization structure functions coming from this part 
of the propagator. Firstly, the description of the spin 0 part of the graviton 
propagator is given. Next, we work out the contributions from the 4-point 
diagram (” Oa”) and the local part of the 3-3 diagram (” Ob”). Working out the 
contribution from the (de Sitter invariant) nonlocal part of the 3-3 diagram 
(”0d”) comprises the largest part of this section. In the end the full results 
for F{x; x') and G(x; x') are given. 

Many of the tensor structure contractions performed in this section were 
checked with the help of the xTensor part of the tensor calculus package xAct 
[99] written for Mathematica. 


4.1 Spin 0 Part of the Graviton Propagator 

The spin 0 part of the graviton propagator depends on the gauge parameter 
b introduced in the gauge condition (|3]). Instead of working with b, we hnd 
it more convenient to use parameter /?, dehned to be 


Db-2 


b-2 ’ 


(134) 


Note that for the range of gauges b > 2 considered in this work, /3 is a 
monotonically decreasing function of b, with (3 > D. 

The spin zero part of the graviton propagator can be written as two 
projectors acting on a scalar structure function [SB] . 


A° 


x; x') = 


-2/32 


■X -x'Pf,^{x') X iAwNN{x]x') 


where the projector P^y(x) is a 2nd order differential operator. 


(135) 


Vfxu{x) = 


9^2v{.x) 

/3 


□ - {f3-D)H^ 


(136) 


The spin 0 scalar structure function iAwp^^^x] x') is obtained by inverting 
twice the kinetic operator for an M| = [f3 — D)H‘^ scalar, and once for an 
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M| = —DH^ scalar. The order of these inversions is irrelevant, we hnd the 
following one convenient, 


C — (/5 — D)H^ x') 

n — {I3 — D)H‘^ iAwNix;x') 
n + DH"^ iAiv{x;x') 


iAwN{x]x') , 

iAw{x]x') , 
i5^{x—x') 


(137) 

(138) 

(139) 


Since one of the masses in the scalar kinetic operators is tachyonic, the scalar 
structure function iAw^Nix] x') must break de Sitter invariance. We sepa¬ 
rate it into a de Sitter invariant and a de Sitter breaking part. 


iAwNNix; x') = WNN{y) + 5WNN{y, u, v) , (140) 


the details of which are given in Appendix C. 

The de Sitter breaking part 6WNN of the scalar structure function (IldOp 
is given in fl248p . The projectors in fll35p . when acted on this de Sitter 
breaking part give zero. 


'P^^,{x) X Vpa{x') X 


6WNN{y, u, v) 


0 , 


(141) 


which means that the spin 0 part of the graviton propagator is de Sitter 
invariant for gauges b > 2 (even though the scalar structure function is not). 
Therefore, we can drop the de Sitter breaking part in fll35p . and write 


A° 


x; x') = 


- 2/32 

{D-l){D-2) "" 


X Vpa{x 


X 


WNN{y) . (142) 


By making use of relations (1256112571) from Appendix C, the action of the 
projectors can be written out as 


i 



0 

per 


(x; x') 


2/?2 

{D-l){D-2) 


DpDMpD'^WNN{y) 


w 

{D-l){D-2) 

2 

{D-l){D-2) 


gp,{x)D'pD'^WN{y) + gp^{x')DpD,WN{y) 
9,^u{x)gp„{x')W{y) . 


( 143 ) 
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i 

[(£>-l)(D-2)] ^ ^t{y) 

1 

2WNN" 

2 

AWNN'" 

3 

WNN"" 

4 

(2 y)WNN''' 2WNN'' 

5 

{2-yfWNN" - {2-y)WNN' - p-y)^ + 


Table 9: The coefficient functions Ci{y) in the expansion fll44p of the spin 0 
graviton propagator in the tensor basis given in Table [3l 


It is convenient to expand this propagator in a de Sitter invariant tensor basis 
dehned in Table |3l 

5 

{x]x') = ^[^,vVp^{x■,x') X Ci{y) , (144) 

i=l 

where the coefficient functions in this expansion are given in Table |9l 




4.2 The 4-Point Diagram 

The contribution to the vacuum polarization from the spin 0 part of the 
4-point diagram (middle one in Fig. [1]) is 


P.JTV 


{X]X') = dM 


AO 

OL^^'^S 


(x; x) 5 { 


^'x — x'] 


(145) 

Here the spin 0 part of the graviton propagator was substituted in the full 
contribution (El), and the 4-point vertex function is given in flT^ . The 
coincidence limit of the spin 0 graviton propagator is 




{x]x) = 4(aFf)V(7^<5)/3C'i(0) + (aTr)V/3^7<5C'5(0) , 


(146) 
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where the coincidence limits of the relevant coefficient functions from Table 
[3 are 


cm 


cm 


T{D+2) 

( 47 r)f "" 4(Z}-l)(D-2)r(f+ 2) 

f r(f)r(i-f)r(i^+Mr(^-M 

\ T{D+2) r(|+Mr(|-M 


X 


1 + 


A 

2b 


'N 


V’ 


D+3 




f D + 3 


■+bN 


+ '0( -^-bN 


2cm + 



T{D + 1) 

(to)?'' (B-i)(£'-2)r(ai) 

£>+i r(f)r(i-f)r(i4i+6A,)r(ai-M 

D r(Z)+i) r(i+6A,)r(i-M 



-CM— 

2bN h V 2 







(147) 


(148) 


They are easily calculated from explicit forms of scalar propagators fl235p 
and fl240p . and expressions fl2511l253p in Appendix C, where, by rules of di¬ 
mensional regularization, all D-dependent powers of y vanish at coincidence. 

The contractions of the tensor structures of the spin 0 graviton propagator 
coincident limit fll46p with the 4-vertex flT^ are 


so that, 

jjfj.uKXay'ySj^ 




{x- x) = 


ay^-yS 


r]^ T] — 7]^ 7] 


x{ [12-(B-4)'^]C,(0) + hLJ(K£_4)c,(0) 


(149) 

(150) 


(151) 
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From here it is straightforward to see that the vacuum polarization is in the 
form given by (ITSD and (122]) in subsection 2.2, 


■‘''‘-4pt0 


{x;x') = a,a;|2p'^[V’"xFo,(a;;a;')+2p^[V’"xGoa(^;:r')(152) 
where the structure functions are 


FQa{x]x') = '^i6^{x — x') 

x|-6hf^Ci(0) 

GQa{.X]x') = 0 . 


, {D-A) . 




H^C,{0) 1^, (153) 
(154) 


In the D = A limit the first term in the brackets in fll53p diverges as 1/(Z1—4), 
the second term is hnite, and the third term vanishes in this limit. This is a 
consequence of the coefficient functions fll47p and fll48p diverging as l/(i2-4). 
Comparing with fl281l29p we see that the entire contribution of fll531ll54p can 
be absorbed into the counterterms by choosing the coefficient 


C-oa = -^<1 -QH^C,{0) + V- (0) 


{D-6){D-A) 


{D-Af 


(155) 

Therefore, the renormalized 4-point contributions to the structure functions 


are, 


Fr = 0 , = 0 . 


(156) 


4.3 Local contributions from the 3-3 diagram 

The ”0b” local part of the 3-3 diagram (left one in Fig. [1]) descends from 
isolating the delta function coming from two derivatives acting on the B- 
type scalar propagator as in (1521) . which we reproduce here (with all indices 
lowered), 

DpD'^y X DxD'^B —)■ -^^^r]p„6l6^i6^{x-x') . (157) 

CL 

and by making the substitution in (I43p . 


IJ.JTV 

J-'-SptOb 


{x]x') = dM jj^u.eaMS- 


ao 

alS^-yS 


(x; x) 6^{x — x') 


(158) 
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Note that this contribution has the same structure as the 4-point diagram 
contribution fll45p . The effective 4-vertex jg constructed by con¬ 

tracting the tensor structure in fll57p with the two 3-vertices in fl43p . 


jjliuK9a0'y5 _ ^/jpreAa/3^ ^O^O-^ucrdtpyS 


Contractions of this effective 4-vertex with the tensor structures in the coin¬ 
cident graviton propagator fll46p are 

, (160) 

= -(/l-4)24^ry"l[-5jl . (161) 

Using the identity 

^ (162) 

we can write the contribution to vacuum polarization as described in subsec¬ 
tion 2 . 2 , 


x;x') = d^d'J 2p'‘[V’''x^ob(2:;a:') + 2 p^[V’''xGofe(x;a:') I, 


flJTl' 

-‘-‘-SptOb 


where the structure functions are 
Fob{x-,x') = K^a^~‘^i6^{x — x') 

x|(3Zl-4)iy^Ci(0) + 
GQb{x]^') = {x — x') 

x|-(3T)-8)i7^C'i(0) 


4^ /r\\ 1 (-^ 'I) 117-4 




{D-Af 


H^C,{0) , 


(163) 

(164) 

(165) 


and C'i(O) and ^ 5 ( 0 ) are given in (1147111481) (recall that they diverge as 1/(T*— 
4) in H —)■ 4 limit). Again, we can completely absorb these contributions 
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into the counterterms 02811291) by choosing the coefficients 

Cob = -^|(3D-4)g^Ci(0) + , (166) 

ACo, = ^|(3D-8)hr2c'i(0) + ^^^hf2^5(0)| , (167) 

making the renormalized contributions to vacuum polarization scalar struc¬ 
ture functions vanish, 

Fr = 0 , GIT = 0 . (168) 

Note that had we not introduced a noninvariant counterterm in fl271) we would 
not have been able to remove the divergence in 0165|) . 


4.4 Nonlocal contributions from the 3-3 diagram 


The non-local ”0d” contribution from the 3-3 diagram derives from making 
the replacement for the two derivatives acting on the S-propagator, 


m expression 


DxD'^B — DxD'^y x B' + Dxy D'^y x B” , (169) 

. This makes the bi-tensor density dehned in fflSj) . 


flKrpuB 

-'Od 


[x;x ) = 


K 


27/2 




X 


D,D '0 y^B' + D,y D’^,y x B "]. (170) 


We use the recipe developed in [6T] and outlined in subsection l2.2l to hnd the 
”0d” contribution to F{x] x') and G{x] x') structure functions. The hrst thing 
is to calculate coefficient functions fi from fl23p . We do this by contracting 
the tensor structure of the graviton propagator with the vertices and the 
tensor structure of the photon part. The contraction rules can be found in 
[95] . This contribution is de Sitter invariant so only fi and /2 can appear (the 
de Sitter breaking part ”0c” is zero in this case), and we hnd it convenient 
to express them as a sum over spin 0 graviton coefficient functions dehned 
in Table E] 


5 

fi{y) = -f^‘^H^^{^^i)i{y)xGi{y) , f2{y) 

i=l 


5 


= X](Av? 2 )i( 2 /) X Gi{y), 

i=l 

(171) 
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i 

(A(pi)i(i/) 

{A^2)i{y) 

1 

-2{5D-13)B' B' 

+2{D-2){2-^y)C' 

2{3D-8)B" + {2-y)B'-3{D-2)C' 

2 

+(D-2)(4imj2)c> 

A{D-3)B' + {4:y-^^)B' 

+ {D-2){2-^)C' 

3 

-liAy-y^fB' 

liiy-y^B” 

4 

-{D-4){4y-y^)B' 

-2{D-A)[{2-y)B'-iD-2)C'] 

5 

-^{D-A^B' 

liD-AfB" 


Table 10: The functions {Aipi)i{y) and {Aip 2 )i{y) defined in equation fll7ip . 
Function B{y) is defined in fl204p . and C'{y) in fl36p . 


where factors are extracted for convenience. The coefficient functions 

(A(^i)j and (A(^ 2 )i are given in Table [Till where identities and have 
been used to simplify the expressions. 

The next step is to construct the master structure functions dehned in 
Tabled] Because of de Sitter invariance T(j/,M,n) has to be zero, and <F(j/) 
is, according to Table HI and expression (I17ip . 


$od(l/) = [^iD-l){Aif^)i-2{2-y){Aif,)[ + D{2-y){Aif2)^ 

i=l ^ 


+(%-|/^)(A</92)- 


X + 


-2i2-y){Acp,), + {4y-y^){A^2)i xC' .(172) 


The results for coefficients of Ci and C*' in this expansion are presented in 
Table dH Plugging in Ci{y) from Table |9] and coefficient functions from Table 
dH into fll72p we get for the master function, 

$od(2/) = X |l6(4|/-|/2)lFiViV"" + 48(2-2/)fTiViV"' 

+ 16(Zl-5)lTiViV" + ^(Zl-4)^| 

+ t(1l-lT "" {-(42/-2/')'WWA""' + 2(Zl-8)(2-y)(4y-|/2)lT^^ 
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i 

2(i>-l)(A¥.i)i-2(2-i;)(A.pi)' 

+0(2-i/) ( A(/32 )i-K4?y-i/2 ) ( A(^2 )' 

-2(2-|/)(A(pi)i + (4|/-|/2)(A(^2)i 

1 

2{2-y)C'-8{D-3)B' 

8{2-y)B'+{Ay-y^)C'-A{3D-A)C' 

2 

4{D-2)(2-^)B' 

+3(4y-i/2)C'-8C' 

iAy-y^)[AB'-{2-y)C'] 

3 

-A{2-y){Ay-y^)C' 

-{Ay-y^yC' 

4 

A{D-A)[{2-y)C'-AB'] 

2{D-A){Ay-y^)C' 

5 

0 

-{D-Ayc 


Table 11: The coefficient functions of Ciiy) and C'(|/) from expression fll72p . 
An overall factor of {D — 2) has been extracted from each term. 


{D^-18D + 70){Ay-y^) - 4:{D^-6D + 32) 


WNN'' 


+ {D-6){3D-U){2-y)WNN" - {D-AfWNN’ 


2 o WN"' 

--iD-A)iAy-y^)^ 


2 WN" 

+ -{D-4){D-6){2-y)^ 


2 WN' {D-4f W'\ 

m 13^ S' 


(173) 


Using the identity for the d’Alembertian operator acting on a scalar function 
that depends only on y, 


□ 

IP 


S{y) 


{4y-y^)S"{y) + D{2-y)S'{y), 


(174) 


and derivatives of it with respect to y, we can rewrite the above expression 
for master function as 


^od{y) = 




[wdl 

^WNN 

1 ^ 



2{D-1) 

+A8{D-1)WNN" + IQ- 


16{D + l){2-y)WNN" 
(D-A) WN" 


/3 




-;-— X < -dy 

2{D-l) ' ^ 




WNN 


+ A{D-2){2-y)d\ 


VH^ 


WNN 
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+A{D-2){D-^)dy 


^WNN 


-^{D-4)dy 


□ PWVI 

IpIp 


-IQ{D-2){D + 1)WNN''' - A{D-2f{D-l){2-y)WNN" 

4 WN" 

+4{D-2f{D-l)WNN' + -{D-4){D-2){2-y)^ 


4 WN' 


{D-4fW' 

^2 


(175) 


Next, using identities (I256il257ll from Appendix C, we can further simplify 
this expression, 


^od(i/) = X <; -i6(zi+i)(2-i/)iyiviv" 


/3 + (2A)-3) 


WNN" + 16 


1 + 


P-4) 


2 (^- 1 ) 

+ 16 

' 2(71-1) 

+4(11-2) \ 13-{D'^-20 + 2)] {2-y)WNN'‘ 


f3 J 772 


WN" 


X -16(71-2)(71 + l)l+AiV" 


- 2 (271^ - 971 +12)/? + (71 - 4) ^ 


1+A7V' 


+4(71-2) 


1 + 


71-4 




( 2 - 2 /) 


WN" 


/?-2 (71^-571+ 8) + 


772 

(71-4)2] WN' 


1 + 


71-4 




=+_±m_ 2 f+ 

Hi + ’ Hi 


(176) 


As argued in subsection 3.5, we can set 71 = 4 in any part of the structure 
functions F{x;x') and G(x;x') that diverges less strongly than l/Ax^, as 
Ax —)■ 0. Since, according to fl241l25p . F is an integral of $ with respect to 
y, and G is a double integral of $ and 4/ with respect to y, in the master 
functions we can set 71 = 4 in any of the parts which diverge less strongly 
than y~^ as // —?■ 0 (recall that y = H'^aa'Ax'^). Also, we can throw away 
terms containing (71 — 4)^. Taking into account the expansion of B' and C 
given in (1206112071) in Appendix A, and functions W, N and NN defined in 
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Appendix C, we get for the master function 


2 TT2D-2-p2( O'- 


^od{y) = 


k^H 






y 


D 


+ 


,D-1 


+ 


y" 


S’ 


where the two D-dependent coefficients are 


IdD) = -- {D-2)f-MD-i)fJ + 


2(C-4)- 


(177) 


(178) 


(B-1) J ’ 

12 (D) = i[/3"(5-/3)-^(4/?2-83/?+180)(B-4)] +o((B-4)2),(179) 

The de Sitter invariant function Af{y) is dehned to be 


3 5 —// 5—// 5 —/// 20—/// 

J ^{ y ) = ~ 3^2 + 32/^2 “ -^^3 -- NN ^ 


,/3(/3-10) 

H-- yNNi — 


.// /3(2/3-5)-_// 5/3 


20/3 
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6 8/3 


r 

f/S 


[10 1 

I 

L8 J 

iVi + 

- 


+ 


^i^-4 

(9 ■ 


NN, - —yNN, + g 

-ITpff 

2/^1 


„ 10 —/// 40—/// 

N 2 - JVN2 + jN2 


NN^ 


(180) 


where the dehnitions in D dimensions of functions Ni{y) and NNi{y) are 
given in fl259il260]l in Appendix C, and the limit H —)■ 4 in fl26ip . which is 
taken here. The power series representation of this function is 

?2 


»r/ 7 r 5 V- n 


7l„lndl +-B. 


(181) 


n=0 


where the coefficients are 

T{^ + bN+n)T{^ — bN+n) 


Qn = 


(182) 


4”+i(n+l)!(n+2)!r(i+6„)r(i-M ’ 

n{n — l)l3‘^ — 4(n —l)(3n + 2)/3 + 40n(n+l) , (183) 


^ (n+ 1 ) 

8 (n + 3)(?T.+4)/3 


Br, = 


(n+ 1 ) 


8/3(n+3)(n+4) . 

X 


n(n —1)/3^ — 4(n —l)(3n + 2)/3 + 40n(n+l) 


'il^[-+bN+n) + 'ilj[--bN+n) - '0(n + 2) - '0(n+3) 
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8(n+3)2(n+4)2/3 


/?2(n^ +14^3+ 3771^-12) 

-4/3(3n^+42n3 + 125n2 + 52n-22) 
+40(n+l)(n^ + 13n^ + 36n+12) . 


Since the master function 4 /( 7 /,M,n) vanishes, according to 
have for the vacuum polarization structure functions, 


Fodix;x') = —2{aa')^ ^xJ 


God{x;x') = {aa') ^xJ 


^^od{y) 

^\D-l)%diy) + y^Uy) 


(184) 
we 

(185) 

(186) 


Performing the integrals above gives the following, 

.2r2fD^ r 2ii{D) 2i2{D)H^ 


Fod{x;x') = 






+ 


K 


(P)-l)aa'Aa;2^-2 {D-2)^x‘^^-^ 

2^6 

^ ^ y-^M{y) 


47r^ 


[aa 


(187) 


God{x]x') = 


^^r 2 (f) 

8ti^ 


ei{D)H^ 


{D-l){D-2)Ax^^-^ 


K' 


iRd 


87r^ 


(aa 


f\2 


2F 


y M{y) +1 y Af{y) } , (188) 


where we have plugged in the dehnition y = aa'H'^Ax^. 

Next we need to localize the divergences using relations fll25p and fll26p . 
which we reproduce here. 


4:71° H6^{x — x') (?2 


^^2D-4 

1 

Ax 2^“2 


2(P)-3)(P)-4)r(f-l) 4 L 

47 rT fd^~^d'^i5^{x—x') 


\n.{n‘^Ax‘^ 


Ax'^ 
d‘^ \\n.{y?‘Ax^) 


4(P)-2)2(P)-3)(Zl-4)r(f-l) 32 


Aa ;2 


(189) 

(190) 


The divergences in (I187p and (11881) are absorbed into the counterterms (l27jl 
by choosing the following coefficients 


= 






n^MD) 


32vrf (/l-l)(Zl-2)(Zl-3)(Zl-4) ’ 


C = {3D-8)G: 


Od 


levrT {D-3){D-4) ' 


(191) 
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Note that, since ^i(-D) = Co{b)/2{D—2), where Co{b) was dehned in ([7]), the 
(7°'^ coefficient indeed coincides with the spin 0 part of the coefficient (I5D]) as 
inferred from the flat space limit [39] . What remains is to perform explicitly 
the integrals in second lines of fll871ll88|l . 


= l\y-V(y) 
0 L 


6 dl3 


<?o^o 
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Qo{^o + Bo) qiAi 2(y 
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(193) 




■A^ciy) = ‘^1 y +i y ^Afiy) 


6 9/5 
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-?o(2v4o+5o)ln(^|j + qiAiy\v?{^ 
+qi{2Bi- Ai)y\n{^ +qi{Ai-Bi)y 
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^ (n+l)(n + 2 ) 


(n+ 3 )y 4„_|_2 In + (n+ 3 )-Bn +2 
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71+2 


(194) 


(n+l)(n+ 2 )‘ 

The following integrals have been used to calculate the two functions above, 

,n+l I- 


2 /" In 
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(n+l) L 4/ (n+1) 


(n^-1) 
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= ll+f^ 
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+ -ln^f^ 

2 V4 
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y 


-1 


= In 


(195) 

(196) 

(197) 

(198) 


where the choice of integration constants is immaterial since any dependence 
on them drops out when the derivatives in fll 8 l) are acted with. Here we 
list explicitly just the hrst few coefficients in fll93til94p of the most singular 
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terms in ?/ —)■ 0 limit, which will be relevant for quantum-correcting Maxwell’s 
equation ([T]), 
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{Qo^o) — 




(/3-5) 
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9 , . . (/^-5) 


(199) 


(/?-2)(9/3-86) 


4/^2 

, (/g-5) 
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(/3-6)(/?-4) 
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, (/5-5) 
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512 3840 
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(/3-6)(/?-4)/-l 
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--■T 


( 200 ) 


( 201 ) 


192 V26jv^ 

where 67 V = [25/4 — Even though the rest of the coefficients can 

be calculated from fll821ll84p . and the remaining series in (11931) and (11941) 
summed into generalized hypergeometric functions, these will give irrelevant 
contributions when used in ([ 1 ]), so we do not do it here. 

The renormalized contribution to the structure functions is 


KTiy) = 


2b- 


1^/ _ 


6-2 7 \487r2' 


aa' 


b-8\ 


/ 6 - 8 \ 

>) 


6-2/ 727r2 
^2 rln(/i2Aa:^)" 


iS^(x—x') 


.6-2/ 5767r4 




3847r"^ aa' I Ax^ 
ln(/i2Ax^)l 


Aa:2 




aa'fMFiy) \ , ( 202 ) 


^ren 

^Od 


iy) = 


2b-l\^( 


6-2 


T2 r 


247r2 
^2772 

~1927r4 


1 — ln(a) i5'^{x — x') 




ln(p2Aa;^)i k^ 77® /.o^r / \ i 

^ ^ +^^(aa')WG(2/) M203) 


Aa;2 


127r4 


Since the renormalized ”0a” and ”0b” contributions (11561) and (11681) vanish, 
the ” 0 d” contribution above constitutes the full contribution to the vacuum 
polarization structure functions from the spin 0 part. Note that there also 
might appear hnite contributions of the form (12811291) . 
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5 Discussion 


We have evaluated the one graviton loop contribution to the photon vacuum 
polarization x') on de Sitter background in the 1-parameter family 

of exact, de Sitter invariant gauges (j3]). The result is represented in terms of 
two structure functions F{x;x') and G{x;x'), whose relation to i[^W]{x;x') 
was dehned in equations fflSj) and fl2^ . Our graviton propagator has a gauge 
independent, de Sitter breaking spin two part and a gauge dependent but 
de Sitter invariant (for b > 2) spin zero part. Each part makes distinct 
contributions to the structure functions. The spin two structure functions 
are fll32p and fll33p : the corresponding spin zero results are fl202p and fl203p . 

The point of this exercise was to check the conjecture [2S] that the leading 
secular dependence of solutions to effective held equations such as ([T]) might 
be gauge independent. The full solutions certainly contain unphysical gauge 
and held variable dependent information because the hat space limit does 
inn]- However, they also contain physical information because one can use 
them to construct the hat space S-matrix |13]. What we need is a hlter to 
distinguish physical ehects from unphysical ones. The S-matrix provides this 
in hat space, but there is as yet no analog for cosmology. There is no ques¬ 
tion that such a hlter exists because astronomers are measuring something. 
Identifying what theoretical quantity represents these measurements is one 
of the central problems of cosmological quantum held theory [FT] . 

Note that the vacuum polarization will play an essential role whatever 
is the outcome. If the conjecture proves to be correct then one can extract 
physical information from the leading secular dependence of solutions to ([T]). 
If the conjecture proves false then one must resort to some form of gauge 
invariant Green’s function, which would inevitably consist of the expectation 
value of the held (which is what solving the ehective held equations gives) 
plus some extra terms to hlter out the gauge dependence HU. 

Unfortunately, we are not yet in a position to answer this fascinating 
question. Had this computation produced the same kinds of terms, with 
diherent numerical coefficients, as occur with the noncovariant gauge [IT] , 
it would not have been necessary to do much work to solve equation ([T|). 
In that case we would simply have read oh the result for each term from 
the previous solutions |15],|16]- Those analyses show that the largest ehects 
derive from the ln(a)z5^(a; — x') part of F(a:; x). It is interesting to note 
that the coefficient of this term is independent of the parameter b, although 
it does not agree with the noncovariant result. We hnd a coefficient of 
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in equation fll32p whereas the coefficient was in equation (136) of the 
noncovariant gauge analysis HU. However, our computation also produced 
some quantitatively different terms at the end of expressions fll32p . fll33p . 
(12021) and fl203p . Each of those terms contains multiplicative factors of {aa')^ 
which should compensate for the measure of the d‘^x' integration. If these 
terms acquire an extra logarithm they can contribute just as strongly as the 
local logarithms. So there seems no alternative to working out another set of 
complicated integrations for these new terms, multiplied by the classical held 
strengths for dynamical photons and for the response to charges and currents. 
Key questions are whether or not the leading secular effects depend on the 
gauge parameter b, and whether or not they agree with the noncovariant 
gauge. It would also be interesting to work out the spin 0 structure functions 
for b < 2. 

One novel feature of our computation is the need for a noninvariant coun¬ 
terterm, despite our use of an invariant regularization and a de Sitter invari¬ 
ant gauge. As was explained in section 2.3, the problem arises from the 
time ordering of the hdAdA interaction and from the fact that the coin¬ 
cident graviton propagator contains a logarithmic divergence (the famous 
“tail term”) which gives rise to a factor of 1/{D — A) in dimensional reg¬ 
ularization. There is nothing we can do about the derivative interactions 
of quantum gravity or the logarithmic divergence of the coincident graviton 
propagator on de Sitter. We might impose a covariant ordering prescription 
on the interactions but this is problematic in view of the need to keep the 
effective held equation ([1]) real and causal by using the Schwinger-Keldysh 
formalism |66l |67l |68l |69l [TOl [TH [72l [73l [71]. 0 The only alternative is to 
resort to the same noninvariant counterterm — the term proportional to AC 
in expression fITTp — that was needed with the noncovariant gauge |41j . 

One major spin-off from our work is the great simplihcation that was made 
in section (3.1) for representing the spin two part of the graviton propagator. 
This should facilitate checking the gauge dependence of other one loop com¬ 
putations involving gravitons such as the fermion self-energy [19115I11521150] 
and the graviton self-energy mm- 

The photon propagator identities (13811391) and their graviton analogs (IMD 
and (l56P are strikingly similar. Equation (l3^ expresses the photon propaga- 

^Note that our wish to use the in-in formalism is only the motivation for employing 
time-ordering. Once that is done, the ultraviolet divergences of the in-in formalism are 
the same as those of the in-out formalism. 
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tor as a longitudinal term plus the tensor —■^dpd'^y{x; x') times the -B-type 
propagator, whereas equation fl56ll gives the (de Sitter invariant part of the) 
graviton propagator as a collection of traces and gradients plus the tensor 
X d'^-^duy times the (invariant part of the) A-type propagator. So 
the tensor structure of the graviton propagator is the square of the photon 
tensor structure, up to terms which drop out when contracted into the ap¬ 
propriate polarizations. This seems like the remarkable insight that on-shell 
gravitational scattering amplitudes are essentially the squares of gauge scat¬ 
tering amplitudes [100111011110211103] . This has usually been thought to arise 
from simplihcations that occur from going on shell but perhaps it is partly 
due to working in exact gauges. 
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5.1 Appendix A: Photon Propagator Functions 

The propagator for minimally coupled scalar with mass M| = [D — 2)if^ 
obeys equation fl35|l and consists of a de Sitter invariant function of y{x;x') 
called iAsix-jX) = B{y). It is closely related — by equation fl36l) — to the 
propagator for a minimally coupled scalar of M| = 2{D — 3)if^, which is 
called iAc{x;x') = C{y). Their expansions are. 


B{y) 


C{y) 


fr(f)/4 

(47r)f lf-lV|/ 
oo - 

n=0 - 

fr(f)/4 

(47r)f If-lU 


r(?7,-i-zi—2) 

r(n+f) 



r(n+f) /yy-f+2' 
r(n+2) Uy 


(204) 
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(?7- + l)r(77- + -D — 3) /|/\”■ (n •^+3)r(n + ^ 

f>+i) U/ r(n+2) Uy 


n-f+2 


n=0 *- 


r(n+2) 


. ( 205 ) 


They obviously agree in D = 4 spacetime dimensions. Our work requires 
only the first two terms in the expansions of their derivatives, 


B'{y) = 
C\y) = 


dvr 2 

A 2. 
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-^+‘^^ + 0{(D-4:)yO 
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y^ 8{D-2)y'^ 
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{D-6){D-4:) 


^ + OUD-A)y^ 
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(207) 


5.2 Appendix B: Spin Two Propagator Functions 

The propagator for a massless, minimally coupled scalar obeys equation 
and consists of a de Sitter invariant function of y{x;x') plus a de Sitter 
breaking logarithm |96], 


iAA{x-,x') = A{y) + k\n{aa') 


The expansion for A{y) is. 
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T(D-l) 
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A{y) = 
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(dvr; 


r(f)74\T-i r(f+i)/4\#- 


D ] D 


-l\y 
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+ Ai 


>^ p(n+T)-l) /yy r(n+f + l) /|/\ 
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where the constant Ai is. 


,(209) 


Ai — 


r(y^-i) 

r(f) 


-V'l 1-^ ) TV' 


D-1 


+ + • 0 ( 1 ) 


( 210 ) 


The full, H-dimensional series expansion for a:') has been 

derived [95] but we here require only the de Sitter breaking part. 
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It is simple to act the scalar d’Alembertian on functions of In(aa'), 


□ /^ln(aa')j =—H^i.{D — l)f'(\n{aa')^ +/" ^ln(aa')j | . (212) 


Because = iAAAABB can use relations (1461147p . (121 ip 

and (I212p to conclude, 


n-{D-2)H‘^ UiA^AABBi^', x') = iAAABix] x') 
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X < - ln^(4aa') + 
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D-l 
D-2J\ 


{D-2){D-l) 
ln(4aa') + constant > (213) 


where the constant is, 

D 
D 

(214) 

This completes the demonstration of equation (ISSp . 

Because our structure functions are at most quadratically divergent we 
require only the leading two terms of J{y) to be kept in D dimensions. Rather 
than acting the derivatives of (12131) on the complicated series expansion of 
iAaaabb^x^x') that has been derived [95] it is simplest to construct Jiy) by 
integrating the differential equation iAaab^x-^x') obeys. 


-1^ ^ 2T)-3 

{D-2f{D-lY ■ 



UiAAAB{x]x') = iAAB{x]x') = 


[iAB{x]x')-iAA{x]x')] 


{D-2)H^ 

From the expansions (12041) and (12091) one hnds, 
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7/ 2 ^ 7/ 2 V / 
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where the coefficients are, 
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, ( 218 ) 

(219) 

( 220 ) 


Expression fl212p defines how to integrate the de Sitter breaking In(aa') 
term on the right hand side of (I215p . The de Sitter invariant analog is, 


□ 9 ( 9 ) = (iV-V^)g"(v) + D(2-y)g'(y) 


( 221 ) 


Equation (12211) suggests that we can solve ng{y) = h{y) as a double integral 
with respect to y, 


^ 9 {y) = Ky) 


yiy) 



' -^[(%-|/^)° ^h{y)\ 
(%-i/2)f 


( 222 ) 


However, relation fl222p is neither tractable in D dimensions, nor even correct. 
The problem with tractability is obvious from the H-dependent powers of 
(%— 1 /^). The problem with validity derives from the impossibility of avoiding 
poles at either 1 / = 0 or 1 / = 4 that would introduce delta functions into 
relation (I22ip which are not present in the desired source function h{y). 

We solve the first problem by extracting the two leading powers of y from 
the solution and then taking H = 4 on the remainder. 
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The remainder Ag{x]x') obeys. 
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We solve the second problem by first setting the lower limit to ?/ = 4 on the 
inner integral of fl222p . which means there are no poles at ?/ = 4. Then the 
poles at y = 0 are cancelled by adding a constant times the D = 4 limit of 
the iA^ix; x'). Because expression fl225p consists of terms proportional to y'^ 
and to ln(|) we only need solutions for these two sources, 


—gi{x-x') = y^ 


gi{x; x') = ~y - ^y^ - ^ In(aa') , (226) 
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(227) 


Here Li 2 (^) is the dilogarithm function, 

U2{z) = - ■ (228) 

Jo t 

The de Sitter breaking factors of In(aa') in expressions (I223p . (I226p and 
fl227p obviously do not belong in J{y), 
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Recall that the four 77-dependent constants Oi and Pi are given in expressions 
(I2171l220p . The functions gi{y) and g 2 {y) are just the de Sitter invariant parts 
of ([226]) and ([2271), 
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5.3 Appendix C: Spin Zero Propagator Functions 

The spin-0 scalar structure function iAwNN^x] x'), introduced in (I135p is 
constructed out of two scalar propagators, iAw{x]x') and iA]^{x]x') [95] , 
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They satisfy propagator equations for a massive minimally coupled scalar, 
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The IT-type scalar propagator has a tachyonic mass M| = —DH^, and 
necessarily breaks de Sitter invariance. We split it into a de Sitter invariant 
and a de Sitter breaking part, 


iAw{x]x') = W{y) + 5W{y,u,v) . (234) 


The de Sitter invariant part is 
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D-2 


(dvr 


r( 


D 


-l\y 


4\f-i 


+ 


r(f+2) 






r(#+ 3 ) 


/4 \ 2 


2(f-3)(f-2)Vy 


-3 /y—2 

+ Wi + W2‘^ 


n=2 


r(n+f + 2)(|) 


y \ n -^+2 


T{n+D){ir 

'^[{n-A-\-2){n-A-\-l)(^n+l)\ n(n-l)r(n+ y) J 


OO 

E 


. (235) 


where 


Wi 

W 2 


r(B+i)fB+i] 
r(f+i)l 2 B y 


r(B+i) 

r(f+i) 




^B + l\ 


(236) 

(237) 


The de Sitter breaking part of the W-type propagator is 


6W(y,u,v) 


kl (T>-l)^e“+ (^^^^)(2-y)n-2cosh(n) I, 


k 


T{B- 1 ) 
(47r)f r(f) 


(238) 

(239) 


In general, the iV-type propagator contains de Sitter breaking parts as well, 
but for the choice of gauge in this work b > 2 {l3 > B), it mass is M| = (/3— 
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D)H'^ > 0. Therefore, it is completely de Sitter invariant, iAp^{x; x') = N{y), 


N{y) 


(47r)T 

r(f)r(i-f) “ rr(|+6^+n)r(|-6^+n) 

r(3-f+ n)(n+i)! 

r(^ + &jv+n)r(^-&jv+n) /y 

r(f + n)n! U 



|, (240) 


where 


Bn — 


(D-l) 



//D + l\ 

Ki> + 2-J 

= V( 2 ) 


(241) 


Note that even though the terms in the series in fl240p cancel in D = 4, the 
whole series is multiplied by a factor diverging as 1/(D —4) giving a finite 
contribution in this limit. 

The spin 0 scalar structure function in fll35p is solved for by inverting 
(I137tll39p . the solution of which is 


iAiYNN{x', x') 


iAisf^ix] x') — iAwN{x; x') 

JIp 


(242) 


where 


iAwN{x]x') = 
iANN{x]x') = 


iAj^{x] x') — iAw{x] x') 


(243) 


pm 

x') = x') . (244) 


The functions in fl2421l244p we can split into de Sitter invariant and de Sitter 
breaking parts. 


iA\YNj^{x; x') 

= WNN{y) + 5WNN{y, u, v) , 

(245) 

iAwN{x; x') 

= WN{y) + 6WN{y,u,v) , 

(246) 

iANN{x]x') 

= NN{y) + SNN{y, u, v) . 

(247) 


The de Sitter breaking parts receive contribution only from the de Sitter 
breaking part of the hh-type propagator fl238p . since the A^-type propagator 
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is de Sitter invariant {SN{y,u,v) = 0), 


SWNN{y, u, v) 

6WN(y, u, v) 
6NN{y, M, v) 

The de Sitter invariant parts are 


SW{y,u,v) 

SW{y,u,v) 

pm 

0 . 


WNN{y) 

WN{y) 

NN{y) 


NN{y) - WN{y) 

Jm 

N{y) - W{y) 

pm 

d N{y) 

dp m ■ 


(248) 

(249) 

(250) 

(251) 

(252) 

(253) 


The de Sitter breaking part in fl234p of the W-type scalar propagator is 
not its homogeneous part, but rather it satishes 


□ 

.m 


+ D 


5W{y,u,v) 




(254) 


which can be calculated by acting with a d’Alembertian on fl238p . Further¬ 
more, 


□ 

.m 


+D 


w{y) 


0 . 


(255) 


Using these relations fl2541l255P together with 0242112441) we can derive useful 
identities for d’Alembertians acting on de Sitter invariant functions 02511 
[2531), 


^WNN(;y) = (P-D)WNN(;y) , (256) 

^tW(j/) = ()J-Z))WW(j/) + t^ + ^, (257) 

^/V/V(j/) = {P-D)NN(y) + !^ . (258) 
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It is convenient to define a dimensionless function which is a part of the 
inhnite sum in the A^-type propagator fl24UI) . 


Niiv) = 


r(f)r(i-f 


r(i+Mr(|-M 5 


T{^ + bN+n)T{^-bN+n) 


r(3-f +n)(n+l)! 


T{^ + bN+n)T{^-bN+n) fy\n 


, (259) 


r(f + n)n! V4 

where the sum starts at n = z, and similarly for the AW-type propagator, 


NN,(v) = , (260) 

The D —)■ 4 limit of function fl259p is 

r(|+&7V+n)r(|-&^+n) ^yy 

^r(i+Mra-M(n+i)!n!Uy' 

~ 2z/>(n+2) + V^(^^+6jv+n) + |, (261) 

where 

No terms arise in this limit from expanding the H-dependence of or /3 in 
fl259p (they all cancel). Therefore, the dehnition of NNi fl260p is still valid 
where now (I26ip is differentiated. 
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